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Abstract 

With the use of pionless effective field theory including dibaryon fields, we study the — )• np 
reaction for the laboratory photon energy El!^^ ranging from threshold to 30 MeV. Our main goal 
is to calculate the neutron polarization Pyi defined as Pyi = ((T+ — ct-)/ {ctj^ + o"_ ) , where cr+ and 
0"- are the differential cross sections for the spin-up and spin-down neutrons, respectively, along 
the axis perpendicular to the reaction plane. We also calculate the total cross section as well as the 
differential cross section (t{9), where is the colatitude angle. Although the results for the total 
and differential cross sections are found to agree reasonably well with the data, the results for Pyi 
show significant discrepancy with the experiment. We comment on this discrepancy. 

PACS numbers: 13.40.-f, 21.45.Bc, 24.70.-fs 
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I. INTRODUCTION 



The induced neutron-spin polarization Py' in the 'yd — np reaction is defined as 
Py/ = (o"_|_ — o"_)/(o"+ + o'_), where 0"+ and cr_ are the differential cross sections for the 
spin-up and spin-down neutron, respectively, along the axis perpendicular to the reaction 
plane. Conspicuous discrepancy between the experimental and theoretical values of Py' is a 



long-standing puzzle in low-energy nuclear physics |2|]. Schiavilla [l| carried out an elabo- 
rate calculation of Py' based on the so-called standard nuclear physics approach (SNPA). In 
SNPA, the nuclear wave functions are generated with the use of high-precision phenomeno- 
logical nucleon-nucleon potentials that accurately reproduce thousands of neutron-proton 
and proton-proton scattering data (for laboratory energies below 350 MeV) along with the 
deuteron properties. The electromagnetic current operators in SNPA are constructed from 
meson-exchange diagrams in such a manner that gauge invariance and the low-energy the- 
orems are satisfied. SNPA has been used to calculate a great many electromagnetic observ- 



ables involving lightest nuclei, and its general quantitative success is well known 
for the 7(i — )■ np reaction, the differential cross sections calculated by Schiavilla 



i. As 

for the 



lab-system photon energy, , up to 30 MeV agree very well with the data. For the spin 
polarization Py', however, there is a large discrepancy between the state-of-the-art SNPA 
calculation jl| and the data. This discrepancy gives a further indication of the seriousness of 
the "Py' puzzle", and one is led to ask whether the problem lies with theory or experiment. 

In order to shed more light on this issue, we study here the 'yd — )■ np reaction in the 
framework of effective field theory (EFT). The application of EFT to nuclear electroweak 
processes, pioneered in Refs. 5|-l7|, has made great progress since then with various specific 
approaches and techniques developed along the way. It is to be emphasized that, insofar as 
all the relevant low-energy constants (LECs) are known, EFT can give model-independent 
results, and that the accuracy of these results can be systematically assessed in virtue 
of the well-defined EFT expansion scheme. In the present work we study the 'jd — )■ np 
reaction in the framework of pionless EFT with dibaryon fields js iOl, which has shown 
good convergence behavior in perturbative expansion for a number of low-energy processes 
in the two-nucleon systems llNl3|. We compare the results of our EFT calculation with the 



experimental data, and also with the theoretical results obtained in SNPA [1]. It is hoped 
that the present study will provide useful information regarding the Py' puzzle. 
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This paper is organized as follows. We describe in Sec. |TT] the basic elements such as 
Lagrangians, the definitions of observables, and the electromagnetic operators. In Sec. Illll we 
enumerate Feynman diagrams that appear up to the order under consideration, and evaluate 
their amplitudes. Sec. HV] explains the relation between the amplitudes and observables. In 
Sec. |V]we show the numerical results and compare them with data as well as with the results 
of the previous theoretical work. Sec. I VI I is dedicated to conclusions. 



II. FORMALISM 



As stated, we work in the framework of pionless EFT with dibaryons (dEFT for short). 



Since we follow the same formalism as in Ref. [1C|, we only give its brief summary here, 
relegating details to Ref. We consider two dibaryons, one in the ^5*0 channel and the 
other in the ^Si channel, and denote them by d"^ and c?*, respectively. A dEFT Lagrangian 



for involving an external vector field is given by 



£dEFT — + £« + £< + ^st 5 



(1) 



where is the standard heavy-nucleon Lagrangian for the one-nucleon sector; Cg (£*) is 
a Lagrangian for d'^ (rf*), while Cgt describes d'^-d^ transition due to an external vector field 
(a photon field). We employ the standard counting rules and calculate the amplitude up to 
next-to leading order (NLO). The Cn relevant to our NLO calculation reads 
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where f ^ is a velocity vector satisfying v"^ = 1, and S'^ is the nucleon spin operator. Here we 
choose = (1,0) and, correspondingly, 25'' = (0,a). = d^ — ^r-V^^ — ^V^ = d^ — iQV^'^^ 
is the covariant derivative involving the isoscalar-vector and isovector-vector fields, Q is an 



electric charge of a nucleon and and v^^^ are the isovector and isoscalar field strength 



fii/ 



tensors, respectively. m]\f is the nucleon mass, while /xy = 4.706 and fis = 0.880 are the 
isovector and isoscalar magnetic moments of the nucleon, respectively. 

The dibaryon Lagrangians, Cs and Ct, and the transition Lagrangian, Cgt, are given by 

1 



-t 



W-V + 
W-V + 



1 



4m 



iV 



[{v ■ Vf - V^] + A, 



U 



Vt 



s\{N^Pi"''^N)+\i.c. 
4(Ar^pf^^)Ar)+h.c. 



(3) 



rriNPd 



'^^^ ieijktltjBk, (4) 



tls3B, + h.c.] , (5) 



T^fi = dfj_ — iCV^^^ is the covariant derivative coupled with the external vector field, where 
C is the electric charge of the dibaryon field in units of the proton charge; C = 2, 1, and 
for the pp-, np- and nra-channel dibaryons, respectively. (A^) is the difference between 
the mass of dg {dt) and 2mN. (yt) specifies the strength of dg-N-N {dt-N-N) coupling. 
Pj;^^"'^ and P/ are projectors onto the ""^^o and ^5*1 two-nucleon states, respectively: 

^r^^ = 4^^2r.a2 , Pf'^^ = ^r,a,a,, (6) 



where and (jj are the isospin and spin operators. B in Eqs. (|3]) and ([5]) is the magnetic 
field, B = V X V'^^^, where V^^^ is the external vector field. Li is a LEG representing the 
strength of a V'^^^-d'^-d^ vertex; L2 is a second LEG parameterizing the strength of a V^^^- 
d^-d^ vertex, and tq are the effective range parameters of the A^A^ interaction for the 
deuteron and spin-singlet channel, respectively. 

The parameters, A^ f and yg^t, in Eqs. ([3]) and (jlj) can be fixed from the scattering length 
and effective range for the ^Sq and ^Si channels. Meanwhile, Li and L2 can be determined 
from the low-energy np — > d'j cross section and the deuteron magnetic moment, respectively. 



Hence there are no unknown parameters in the Lagrangian, £dEFT, in Eq- ([I]); see Ref. 10 1 
for further details. 



III. TRANSITION AMPLITUDES 

The Feynman diagrams contributing to our NLO calculation are depicted in Fig. [H The 
transition amplitude A, in the center-of-mass (cm.) frame, corresponding to diagrams (a), 
(b), and (c) in Fig. [T]may be written as 

A = Xi^(^2r2X2^ ■ {[e(d) x (A; x e^^))]XMs + ■ pYes} 

+Xi(^2r?,T2xV i^d) ■ {k X e(^))XMv 

+Xi^(^2r-iT2xV ■ {e(d)e(7) ■ PXev + [e(d) x (A; x e(^))] Ymv} 

+xU2r2xV i^d) ■ (k X e(^)) Yms , (7) 
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(a) (b) (c) 



FIG. 1: Diagrams for — t- np reaction: a double line with a filled circle stands for a dressed 
dibaryon field, a single line for a nucleon, and a wavy line for a photon field. The photon-nucleon- 
nucleon vertex with "x " is of NLO, and the photon- dibaryon- dibaryon vertex with "x " is propor- 
tional to the LEG, Li or L2. 



where ej;^) and q^) are spin polarization vectors for the incoming deuteron and photon, 
respectively, while Xi X2 the spinors of the outgoing nucleons. k is the momentum 
of an incoming photon (which is taken to be in the 2;-direction), and pis the relative three- 
momentum of the two nucleons in the final state, and we have introduced k = k/\k\ and 
p = p/\p\. The coefficients of the terms in Eq. ([7]) are given as 
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where p = \p\, and u is the incoming photon energy in the cm. frame. 



(12) 



IV. DIFFERENTIAL CROSS SECTION AND NEUTRON SPIN POLARIZATION 



We calculate the differential cross section and neutron spin polarization Py' in the cm. 
frame. ^ The differential cross section is given as 

da a pEi ■^1^,2 
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The symmetry factor 5" is equal to 2 in the present case. In Eq. ( IT3l) . a is the fine structure 
constant, Ei = Jm?^ + is the energy of an outgoing nucleon in the cm. frame, and 



P 



(15) 



where is the mass of the deuteron. The total cross section is obtained by integrating 
Eq. (|T3l) over the direction of p. 

To calculate the neutron spin polarization, we introduce the spin-isospin projection op- 
erator. 



(16) 



where n is the neutron spin polarization axis. We follow the convention for coordi- 
nates in [1^, from which we have k = (0, 0, 1), h = y' = (— sin0, cos0, 0), and 
p = (sin 6' cos 0, sin 6' sine/), cos^), where 6 and are the colatitude and azimuthal angles 
in the cm. frame. Inserting the projection operator in the spin-isospin summation of the 
squared amplitude, we obtain 



MS, 



spin 



+2 IXmvI + I^msI — {Xl^.iyYMs + Yl^gXMv) 
+3 [l-{k- pf] [\Xev? + \Yes? - {X^Yes + Y*esXev) 
±ih -{kxp) [{XlyXMv - Xl,yXEv) + {Y*sYms - YIjsYes) 
~ (YesXmv - Xl^yYEs) - {Y^,jyYMs - Y^sYmv)] ■ (17) 



^ If comparison with the experimental data necessitates it, we shall convert them into laboratory-frame 
quantities. Numerically, this conversion is not important in the present case. 
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FIG. 2: Total cross section, a, for the jd — t- np process as a function of the lab-frame photon 
energy El!^^ (MeV). The long-dashed and short-dashed lines show the Ml and El contributions, 
respectively; the solid line gives the sum of the El and Ml contributions. The points with error 
bars represent experimental data. (The experimental data are obtained from National Nuclear Data 



Center (NNDC) web-page 



Noting that, whereas Xmv and Xms are complex, Xev, Ymv, Yms and Yes are real, we 
arrive at a final form for the polarization Py/ as 

= -2 sin e{XEv - Yes) I^Xmy/ {4 {\Xms\^ + \Ymv\'' - '2Ymv ^^Xms) 

+2 {\Xmv? + \Yms? - 2Yms'^^Xmv) + 3(1 - cos^ ^) (\Xev? + \Yes? - 2XevYes)] ■ 

(18) 

Since £' = (cos ^ cos 0, cos 6' sin 0, — sin 6*) and z' = (sin 6* cos 0, sin ^ sin 0, cos 6'), one can 
easily verify that Px' and Pz> vanish in the chosen coordinate system. 



V. RESULTS AND DISCUSSION 

Fig. [2] shows the total cross section, cr, for the 'jd np reaction from threshold to 
Ei^'' = 30 MeV. It is seen that, at low energies, there are a few data points that are off the 
calculated a curve. We remark, however, that the error bars in the figure only represent 
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FIG. 3: Differential cross section 100 a{9)/a at El""^ = 19.8 MeV, where 9 is the scattering angle 



in the cm. frame. The experimental data labeled "De PascaleSS" are taken from Ref. il a /. 



statistical errors, and that the error bars are hkely to become significantly larger when 
systematic errors are included. The data at higher energies, El!^^ > 15 MeV exhibit some 
scatter, but their overall behavior is consistent with the calculated cross sections. Thus we 
conclude that our dEFT calculation up to NLO, which contains no adjustable parameters 
after the two LECs (Li and L2) have been fixed, can reproduce the total cross section data 
reasonably well. 

In Fig. [3] we plot the differential cross section cr{9) at E':^^ = 19.8 MeV, where 9 is the 
scattering angle in the cm. frame; the figure also shows the data from 16|. It can be seen 
that the calculated differential cross section is consistent with the measurement. 

Fig. m shows the angular distribution of the polarization Py/ calculated for E^^'^ = 2.75 
MeV, along with the experimental data. There are two sets of data available in the literature. 
One set (referred to as John61) is due to John et al. [171], and the other set (referred to as 
Jewell65) is due to Jewell et al. ISj. We note that John61 has significantly larger error 
bars than Jewell65. Fig. S] indicates that our results agree with John61 within the large 
error bars except at cos0 ^ —0.75. However, compared with Jewell65, the theoretical curve 
clearly lies below the experimental values for the entire angular range. In fact, this pattern 
of discrepancy between theory and experiment was already discussed in Ref. {l8|, where 
the authors used the theoretical values of Pyi that would turn out to be close to what we 
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FIG. 4: Pyi at E}^^ = 2.75 Me V plotted as a function of cos 6. The experimental data labeled 



"John61 " are from 



j, and those labeled "JewelWS" from UaJ - 



have obtained here. It is to be added that Py' for E'^ 



lab 



2.75 MeV calculated in SNPA |l| 
agrees with our results. Thus the data set Jewell65, which has much smaller error bars than 
the earlier set John61, disagrees with both the latest SNPA and dEFT calculations.^ The 
persistence of the Py/ puzzle suggests the desirability of a new measurement of Py>. 

Fig. [5] shows, as a function of -Ei"*, the Py' for the lab-frame scattering angle 9iab = 45° 



calculated in dEFT; also shown are the experimental data taken from [19|. The figure 
includes the results of the previous SNPA calculation jl| as well; the values labeled "lA" 
correspond to the impulse approximation (lA) while those labeled "MEC" include the meson 
exchange currents. We can see from Fig. [5] that the results of our dEFT calculation com- 
pletely disagree with the data. The figure also indicates that the present dEFT calculation 
gives values of Pyi significantly different from those obtained in the SNPA calculation [l| 
(although, for E!^^ < 10 MeV, the dEFT curve happens to be rather close to the lA values 



in fl. 



Whereas the SNPA values increase almost linearly as functions of E^ (with the lA 
curve increasing more rapidly than the MEC curve), Py/ obtained in dEFT almost tapers 
off around E'i"* 10 MeV. This latter feature worsens the discrepancy between experiment 



We remark en passant that, if we multiply the calculated values of Py> with a factor of about 0.7, the 
scaled results agree with Jewell65. 
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FIG. 5: Pyr at Oiab = 45° plotted as a function of , where 9iab is the scattering angle in the 



laboratory frame. The data labeled "Nath72" are from ]19i ]. Filled triangles and open circles are 
results of I A and MEC in SNPA, respectively /J/. 



and theory, which was aheady conspicuous with the use of the SNPA values of Pyi. 

Fig. Oshows, as a function of -E';^"'', the Pyt at Oiai^ 90_° calculated in dEFT, together with 



three sets of experimental data taken from Refs. |19l42l|. Comparison between theory and 
experiment is hampered by the fact that the data points show pronounced scatter, exhibiting 
at some places even inconsistencies among the data sets from the different sources. It is 
also curious that, for E^^^ > 10 MeV, the data points show rather conspicuous oscillatory 
behavior. We may summarize the current situation with the statement that the average 
behavior of the experimental values of Py'{Oiab= 90°) agrees with the results of our dEFT 
calculation within the very large experimental uncertainties. Here again, it seems desirable 
to have new measurements of Py'. Fig. |6]also gives Py> obtained in a SNPA calculation 1|. 
It is seen that, within SNPA, the lA calculation always gives larger values of Py' than the 
MEC case. The curve corresponding to our dEFT calculation lies between the lA and MEC 
values for E^^^ < 10 MeV, but it approaches the lA results as E^^'' increases. 

In Fig. El we plot the dEFT values of Py> at 9iab = 135°, as a function of_E^''^ For 
comparison, we also show the SNPA values of Py/ for the lA and MEC cases [l]. It is to 
be seen that, for Ei^^ < 5 MeV, the dEFT and SNPA results are close to each other, but 
that qualitative difference appears for El^'' > 8 MeV; the dEFT curve shows slow, monotonic 
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FIG. 6: P,/ at Oiab = 90° plotted as a function of E^^^ . The data labeled "Holt83" and "Drooks" 
are from '2G] and \21j, respectively. Filled triangles and open circles are results of I A and MEC in 
SNPA, respectively Jj/. 
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Filled triangles and open circles are results of lA and MEC in SNPA, respectively Jj/. 



increase, whereas the SNPA results (both lA and MEC cases) start decreasing around E}^^ > 
10 MeV. We also remark that, in our dEFT calculation, the dominant contributions to Py' are 
proportional to sin^ or cos^ 9, which implies Pyi{6iab= 135°) ~ Py'{Oiab= 45°). Fig. [TJ which 
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also includes Py'{9iab= 45°), indicates that Py'{9iab= 135°) ~ Pyi{Oiab= 45°) holds rather 
well. By contrast, the calculation in |lj] does not share this feature, another qualitative 
difference between our results and those in jl|. As regards comparison with experiment, 
because we were unable to retrieve the relevant data from the literature, we cannot make 
direct comparison of the theoretical Py>{6iab= 135°) with experiment. However, to the extent 
that Pyi{6iab= 135°) ~ Py'{6iab= 45°) holds, we can compare our theoretical curve with the 
data for Pyi{6iab= 45°) shown in Fig. [5] . 

VI. CONCLUSIONS 

We have applied the pionless-EFT-with-dibaryon (dEFT) formalism to the '-yd — )■ np 
reaction for the incident photon energy up to 30 MeV. As far as the total cross section and 
the differential cross section are concerned, the results of our dEFT calculation agree with 
those of the latest SNPA (standard nuclear physics approach) calculation by Schiavilla [l|. 
These theoretical values are in reasonable agreement with the data, which at present have 
appreciable uncertainties. On the other hand, for the neutron polarization Pyi , the results of 
our dEFT calculation are found to be significantly different from those obtained in SNPA jl|, 
indicating the sensitivity of polarization observables to the theoretical frameworks used. It 
is noteworthy that, even if we interpret the difference between the EFT and SNPA results as 
a rough measure of the existing theoretical uncertainties, the "-P^,'" puzzle, i.e., discrepancy 
between the theoretical and experimental values of "-Py'", still persists. Our results indicate 
that Pyi obtained in dEFT can exhibit even larger discrepancy from the data than the SNPA 
calculation does, for certain ranges of the scattering angle. 

We remark that, at energies larger than E}^^ ~ 15 MeV, the contributions of final-state 
partial waves higher than considered in our present calculation may become significant, 
and that the SNPA calculation [l] includes these higher partial waves. This may explain 
part of the differences between the dEFT and SNPA results for Pyi. Higher order effects 
within dEFT also need to be examined despite the good convergence property of dEFT 
found previously for many observables. It is to be noted, however, that Christlmeier and 
22( 1 have carried out an N^LO calculation in dEFT for the longitudinal and 



Grief^hammer 



transverse response functions for the (i(e, e') reaction. According to this work, it is highly 
unlikely that the large discrepancy between theory and experiment found for some of these 
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response functions can be ascribed to higher order terms in the dEFT expansion. A similar 
conclusion may hold for P^/, and a calculation going up to N^LO for the 7(i — > np reaction 
seems warranted. We also remark that, even at energies below 10 MeV, the inclusion of 
higher order corrections is desirable in that it will reduce theoretical uncertainties and help 
sharpen the issue of the discrepancy between dEFT and SNPA at low energies. 

We make here a brief comment on the treatment of the internal structure of the deuteron 
in dEFT. The introduction of the elementary dibaryons, ds and rff, in dEFT might give the 
impression that the deuteron structure has no place in dEFT. It is to be noted, however, that 
a photon coupled to the intermediate nucleon (diagram (b) in Fig. 1) gives rise to momentum 
dependence in the deuteron form factor, and thus the structure effects subsumed in the form 
factor can be accommodated in dEFT. In Ref. [10], the electromagnetic form factors for the 
deuteron were calculated in dEFT up to N^LO, and the differential cross sections for e-d 
elastic scattering were computed with the use of these form factors. The functions, A{q) 
and B{q)^ that represent the momentum dependence of the cross section (see Ref. jlO| for 
the definitions of these functions) were compared with the experimental data and also with 
the results of other theories, and good agreement in the low-momentum transfer region was 
reported. Given the generally good convergence properties of dEFT, we may expect that 
our present NLO calculation incorporates the bulk of the deuteron structure effects, even 
though a possibility does exist that Pyi is a "delicate" quantity that is exceptionally sensitive 
to higher-order terms. In this context also, an extension of the present work to higher chiral 
orders seems of importance. 

It is worth emphasizing that the accurate understanding of polarization observables is 
also important in connection with parity-violating observables in nuclear electromagnetic 



processes 



231]. In the process 70? — )■ np, for example, the neutron polarizations along the x' 



and z' directions vanish with the parity-conserving interactions as mentioned before, but they 
can be non-vanishing with the parity-violating interactions. Theoretical prediction on these 
parity-violation observables requires high accuracies in both the strong and electromagnetic 
amplitudes. A polarization observable which is sensitive to the interference between the 
strong and the electromagnetic amplitudes can be a good testing ground for the reliability 
of parity-violation calculations as well. 

To summarize, our study points to the necessity of further studies, both experimental 
and theoretical, of the spin observables in the 7(i — )■ np reaction. 
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